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ABSTRACT 


“1. Introduction 


Maize (Zea mays L.), the "queen of cereals", is the most widely used cereal in the world's diet and also in the 
manufacture of animal feeds, particularly in sub-Saharan Africa; this makes it the most important cereal in terms of 
production and consumption ahead of rice and millet in West and Central Africa [1],[2]. It is generally preferred to 
other cereals because it has a more uniform nutritional value and provides better structure to the food, and above all 


because of its availability [3],[4]. 


Around two-thirds of the maize grain produced in Cameroon is used in feed mills for the production of animal feed 
[2]. As a general rule, the first stage in the production of animal feed from maize grain is grinding, which is carried 
out thanks to the development of mills and grinders equipped with electric or thermal motors. These mills, which 
used to be imported, are now, for most of the mills used in Cameroon, locally manufactured [5]. These mills or 
grinders are generally hammer mills, and are the most widely used for dry grinding of maize due to their simple 
design, easy operation and low maintenance costs [6],[7],[8]. In these mills, the product to be ground is introduced 
into the grinding chamber by gravity and exits the cylinder when they are small enough to pass through a screen at 
the bottom [9]. Although functional and increasing in demand (due to the cost of purchase), this locally-designed 
equipment does not always meet users' expectations, and does not preserve the organoleptic qualities of the flour 


obtained. 


The process of grinding cereal grains in hammer mills is complex and not fully understood [10]. Currently in 
particle milling, producing a ground product to the desired size with minimal energy consumption is the most 


challenging task for a process engineer [11]. It is important to understand the factors that affect performance 
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parameters such as energy consumption per unit mass and the particle size distribution of the product leaving the 


mill [12],[13]. 


Several studies have been carried out to contribute to this scientific challenge. Some studies are mainly focused on 
investigating the efficiency of the hammer mill in terms of particle size distribution of the output product, energy 
consumption and yield [7],[14] by adopting an empirical approach that consisted in modeling these performance 


criteria as a function of a few defined variables. 


On the other hand, other studies are limited to investigating the kinetics of particle breakage in a small hammer mill 
[15],[16] and the strength characteristics of particle under impact fracture [17]. However, studies evaluating the 
effect of hammer mill parameters on performance criteria are rarely reported in the literature [18],[19]. Note also 
that semi-empirical energy models by Rittinger, Kick, or Bond, have been widely applied to model the energy 
consumed during grinding as a function of product size at the mill inlet and outlet, during grinding in hammer mills 
[14],[20]. Other semi-empirical models used to smooth the particle size distribution curves of the grinding result 
using simplified laws such as the Gaudin-Schuman, Rosim-Rammler-Bennet or law [21],[22]. However, these 


approaches are incomplete, and lack appropriate physical understanding. 


In order to comprehensively predict the complete particle size distribution for a given set of operating conditions, 
population balance models have frequently been used. It is most often used to mathematically describe the crushing 
process [12],[23],[24],[25]. In the following, we present the population balance approach as a tool for mechanistic 


modeling of grinding in hammer mills. 
“= 2, Mathematical Modeling of Hammer Mill 
2.1. Population Balance Model (PBM) 


In the case of grinding, in addition to their ability to predict the particle size distribution of the process, mechanistic 
models contribute to the understanding of the physical phenomena involved. They are based on the fundamental 
laws governing the behavior of the process under study. According to [26], grinding can be considered as a 
succession of unitary stages, essentially consisting of the selection of a material fraction in the first stage, followed 
by the grinding of this material fraction in the second stage, to which is added the classification of this material 


fraction in the case of screen mills. 
The general form of the population balance equation for grinding is described as follows [27],[28]: 


dm(x,t) 


=| Input | -| Output |+| Birth|-| Disappearance (1) 
dt 


In equation (1), m(x, t) represents the mass fraction of particles of size x present in the mill during the time interval 
dt. Birth and Disappearance respectively represent the increase or decrease in particle mass fraction due to the 
grinding process. Input represents particles entering the mechanism, and output represents particles leaving the 


system. 


Based on this kinetic approach, [29] proposed equation (2), representing the discrete form of equation (1), to model 


hammer mills with a built-in sieve [30]: 
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dm,(x,t) ° ° i-l 
a = Min, (Xt) — Mout,i Gt) + Y SiC fo)m(O-SGIm (2) 


where the output is, 
Mout,i (x, t) a C; (x)m,; (t) (3) 


In equation (2), one, (x,t) and jis (x,t) represents on the time t, the particles outflow and inflow rates in the mill, the 
particle was size or fraction i. §,(x) (s") is the rate of disappearance of particles of size i that are ground at time t, or 
selection function. b;;(x) is the size distribution of particles of size j that have ground to size i at time t, or grinding 
distribution function. C(x) is the rate of the mass fraction of particles of size i leaving the mill, and t the residence 


time of the particles inside the mill. 


This equation implies that the quantity of material of size i evolves due to three phenomena: 
e It decreases as size i decreases, expressed by the selection function S,(x). 


e It increases due to the addition of material of size i, resulting from the grinding of larger particles. The grinding 


distribution function refers to the material of size j that fragments for size 1. 


e It also decreases due to particles of size i leaving the mill, reflected by the classification function C;. 


In order to develop a mill model, it is generally assumed that the material passing the screen has remained in the mill 


for an extremely short time [31]. 


This model is based on the determination of these three functions. In practice, this means the evolution of the 
particle size distributions of particles of any size over time. The challenge then is to select an appropriate selection 
function §;, milling distribution functions b; and classification matrix C;, and to make these choices physically 
relevant [31]. There has been some progress in basing these functions on material properties and equipment design 
variables [32],[33], but fracture modeling still remains difficult due to the difficulty of measuring their intrinsic 
mechanical properties, which are linked not only to the fracture behavior of the material but also to the operating 


conditions of the equipment. 


The vast majority of research work that has focused on modeling grinding with a classification system external to 
the mill [11],[14],[34]. Reference [35], [36] and other researchers, build on this concept to develop models for 
simulating and predicting grinding processes. They defined grinding selection and distribution functions, 


continuous in time and discrete with respect to the size classes of ground particles. 


Still on the subject of applications, [10] developed a steady-state population balance model of a laboratory-scale 
hammer mill to predict the particle size distribution for two biomass resources, which are characterized by 
non-brittle fibrous behavior and non-spherical particle shape, analyzing the effects of material properties and 
operating variables. Long before him, [23] explored the modeling of hammer milling of corn, in particular for 


continuous and regular milling in a laboratory-scale apparatus. 
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Reference [9] proposed the dynamic modeling of corn kernel grinding in hammer mills using a population 


equilibrium approach. Corn grinding was tested in pilot and industrial continuous hammer mills, under different 
operating conditions. The results of the pilot experiments were used to adjust the parameters of the selection and 
classification functions, taking into account the effects of rotor speed and screen aperture size. The model was 


validated with steady-state data at both pilot and industrial scale, which also included variations in feed rates. 
2.2. Matrix Form for Hammer Mill Modelling 
2.2.1. Matrix form of PBM 


The equation of the balance approach can be solved numerically in several ways, including using the matrix method 
[37],[38]. The matrix form of the population balance model for predicting hammer mill performance is inspired by 


the structure model adapted from [10]. 


Figure | below shows the hammer mill structure model used for the matrix form of the population balance equation. 


For this case study, the sieve or classification system is incorporated into the mill. 


Dre - impacted Fragmentation : . 
Input Pa ee Impacted Classification Output 
Material h Particles 
m = l 
C— LU 


Recirculating 
Particles r 


Figure 1. Representation of a hammer mill with an adapted sieve by [10] 
In this fig. 1, » is the size distribution of the input material, ; the size distribution of the material after impact and 


before the internal sieve, P is the size distribution of the output product leaving the mill throughout the internal 


sieve and r the recirculating particles that were unable to pass through the sieve according to their size. 


The matrix script adopted for grinding in the hammer mill is therefore: 
p=C-X-[I-U-C)-x] +m (4) 


where I is the identity matrix. C is the classification matrix, with is diagonal matrix whose elements represent the 


classification efficiency for each size interval of material. X is the fragmentation matrix, which summarizes the 
fundamental mechanisms governing particle breakup in the mill, and particle-particle and particle-screen collisions 


are considered negligible. It is presented in equation (5) below: 
X =b-S+I-S (5) 
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whereb is a lower triangular matrix that represents the breakage distribution matrix and S' is a diagonal matrix 


which means breakage rate matrix. 

These relationships are widely accepted by many authors [10],[38] |39],[40]. 
2.2.2. Expression of different function 

2.2.2.1. Selection function 


The values of the breakage rate function also called selection function S$; describe the grinding rate during the unit of 
time. It is strongly influenced by grinding conditions (especially rotational speed) and material properties [9],[33]. 


Its expression is defined by equation (6) [41],[42]: 


Ay 
S,=AN| (6) 
d 


where x,; is the representative particle size in class i. N is the rotational speed of the hammer mill rotor in tr/mn and 
dyo is size reference equal to the maximum corn kernel size. The parameters 4; and A; are the parameters of the 


function to fitted using experimental data. 
2.2.2.2. Breakage distribution function 


The breakage distribution function bj describes the new particle population generated from grinding the parent 
particle. The grinding function is characteristic of the product to be ground. The form for the cumulative breakage 


distribution function proposed by [43], and intensively used by various authors [6],[9],[14], is defined in equation 


(7): 
*pi-a . *pi-a id 
By = o (2) +(1-) (2) aie (7) 
1 i<j 
And the breakage distribution function obtains from the cumulative function is: 


b, = B, — Bij (8) 


where x,;.,1S the representative particle size on the class i-1, x,; is the representative particle size of input material on 
the class j ; g, y et # are the parameters characterizing the breakage distribution function or its cumulative form that 


must be adjusted on the basis of experimental data. 
2.2.2.3. Classification function 


When the particles undergo grinding, they fall onto the sieve, where they are either retained to be ground into fine 
particles, or leave the grinding chamber, depending on their size. The population balance approach defines a 


probability function called the classification function, which translates this step. 


Following [31], [44] proposes a classification function (defined by the C; matrix) defined by equation (9) below: 


6 aie( 2.) (9) 
wad, 
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with 


u=f(N,d,)=a,+a,N +a,d, 


where d, is the size of the sieve mesh, w is the critical size factor that takes into account the rotational speed of the 
hammer mill rotor and the size of the sieve mesh, @, is the adjustment parameter that determines the sensitivity of 
particle size on the efficiency of the classification function. yd, is the maximum particle size that can exit the mill. 
As particle size approaches this maximum size of particles that can exit the mill, the probability of classification 


becomes zero, i.e if x,; > yrds, then C; = 0 (particles are held back to undergo possible fragmentation). 


Two methods exist for determining $;, b and C;. The most widely used, which uses both direct measurement and 
back-calculation from experimental grinding data [45],[46],[47]. And the other so-called "direct measurement" 
method requires tedious breaking experiments on many one-dimensional feeds to determine the parameters without 


resorting to complex non-linear optimization methods. 


Particle size distribution is one of the main sources of information for milling, as it has a considerable influence on 
other processes after milling. One of the aspects to be taken into account when modeling flour particle size 
distribution is the variation in the production of fine particles (sizes below 0.25 mm). The greater or lesser presence 
of these particles is a source of multiple problems, including explosion risks, health hazards and product loss in the 
form of emitted dust [48]. Yet local SME workers are more exposed to these risks due to the use of ageing 
equipment, an inadequate working environment, a lack of awareness of the health risks due to exposure to airborne 


particles [49],[50]. 


Thus, process optimization, based on accurate prediction of particle size distribution at the mill outlet, should aim to 
reduce the rate of fine particle production and cut energy consumption by avoiding over-grinding. And today's 
requirements for the manufacture of processing equipment make it essential to develop research into cereal milling, 
with a view to ensuring food quality that meets the physiological needs of consumers. The most important research 
activity in cereal milling is therefore to obtain a final product with a uniform granulometric composition, a desired 


particle size and minimal energy costs. The aim of our work is therefore to: 
e model the particle size distribution of corn flours obtained from a locally manufactured mill; 


e determine the operating conditions that will reduce the quantity of fine particles produced, while consuming less 


energy during grinding. 

“3, Materials and Methods 

3.1. Material 

3.1.1. Locally fabricated hammer mill 


A pilot hammer mill is manufactured in a workshop in Ngaoundere by a local manufacturer. It is fitted with a rotor 
carrying 16 hammers and a sieve. The diameter of the hammer tip is 470 mm and the width of the crusher is 70 mm. 
The hammers are rectangular, 43 mm wide and 150 mm long. The specifications of the hammer mill used for the 


experiment are listed in Table 1. 
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Table 1. Features of the locally manufactured hammer mill 


Fabrication Locally manufactured 
Capacity 1000 kg/h 
Number per line 4 
Line 4 
Hammer Shape Plate 
Dimensions (mm) _ |4x 50x 120 
Weight (kg) 250 g 
Sieve mesh (mm) 2mm et 8 mm 
Origin Electrical 
Power 7.5 kW 
Energy 
Transmission Pulleys - belt 
Motor speed 3000 tr/mn 
Design Circular 
Grinding 
Width (mm) 250 mm 
chamber 
Inner diameter 600 mm 
Space between hammer and sieve Between 3 to 4mm 
Diameter of the rotor (mm) 390 hammer edge to edge 


3.1.2. Measurement and laboratory equipment 


a K-tronic vibrating sieve shaker equipped with 6 sieves with mesh sizes of 0.25 mm, 0.425 mm, 0.60 mm, 0.85 
mm, 1.00 mm and 1.6 mm respectively, plus a collection bottom. This will enable particle size analysis of the flour 
fractions obtained from grinding corn kernels. A high precision digital balance was used to weigh the maize kernels 


and for particle size analysis. And we use a wattmeter for measuring the power consuming by the hammer. 
3.1.3. Crops material 


White maize, composite variety 'Shaba' (see in fig. 2.) was selected with a moisture content of 13 + 1%. Impurities 
such as plastic strips, stones, maize straw and metal were removed from the maize batch by sieving through a 6.3 


mm wire mesh combined with natural airflow ventilation, and also removed manually. 
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Figure 2. White maize, composite variety “Shaba” 


The maize kernels are collected on the local market. Maize kernels of this variety are generally between 6 mm and 


12 mm in size [51]. Using 6 mm and 10 mm sieves, maize kernels are grouped into three classes, as shown in the 


table below: 
Table 2. Different size classes of maize kernels 
Classes I II I 
size (mm) < 6.3 mm [6.3—10] > 10mm 
3.2. Method 


3.2.1. Experimental procedure 


To develop a model for predicting product size distribution in locally manufactured hammer mills, three categories 


of experiments (experiments A, B and C) were carried out. 


e Experiment A: amass m=4 kg of corn grain of each class is introduced into the mill, with the outlet blocked by 
an unperforated metal plate, for 5 minutes. The variables defined for this experiment are feed rate (Partially Open 


and Totally Open) and rotor speed (2000 and 3000 rpm). 


e Experiment B: operating conditions are the same as for experiment A, but the metal plate is removed and 


replaced by sieves with aperture sizes of 2 mm, 5 mm and 8 mm. 


These two experiments will be used to determine the parameters of the selection, grinding and classification 


matrices. 


e Experiment C: Corn kernels are fed into the grinder without taking into account the different classes. The aim of 
this experiment is to validate the performance model defined by matrices S, B and C. the composition of the 


samples used in this experiment according to grain size is given in the table 3 below: 
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Table 3. Composition of the samples 


Kernel size <6.3 mm [6.3 — 10] > 10mm 
Mass fraction (%) 0.22 0.51 0.27 

Xpi (mm) 4.35 7.94 11.18 

d; (mm) 8.025 


In this experiment, we measured the power consumed. And using (equation 10), we calculated the specific energy 


consumption by the hammer mill. 


E,.. = AP xt 
, M 


(10) 


In (equation 10), P is the difference between the power consumed at no load and the load power consumed, t is the 


time during the grinding has done and M is the grinding maize kernel. 
In terms of operational parameters, three variables are considered: 


e The size of the sieve in the mill (ds), which impacts the product size distribution. Three mesh sizes were chosen: 


2mm, 5 mm and 8 mm. 
e Rotor speed (N), which affects impact energy: 2000 tr/mn, 2500 tr/mn and 3000 tr/mn. 


e Mill feed rate, which influences grain-to-grain interaction inside the grinding chamber. Depending on the habits 
of local operators, two feed rates are chosen, in relation with size of opening hopper feeding: either totally open 


(TO) or partially open (PO). 


All the combinations of these variables give us a total of 20 tests for experiment A, 30 tests for experiment B and 18 


tests for experiment C. 


Particle size analysis is performed by sieving 100 g of ground product using a K-tronic sieve shaker with six sieves 


of respective sizes of 1.6 mm, 1 mm, 0.8 mm, 0.63 mm, 0.45 mm, 0.25 mm and pan. 
3.2.2. Experimental values of S and B matrix 
Experimental values of § and B are determined using the relationships below. 


S values are assumed to be independent of time [52]. The values of the Si matrix are calculated by (equation 11): 
; S.xt 
og| a | 5 (11) 


where 71,5; 1s the mass of unground material (g) in class 1, mg is the initial mass of maize introduced into the mill and 


t is the the estimated grinding time for the mass of 4000 gr. 


For values of bj, corresponding to the fraction of mass of size j that appears in size interval i after primary grinding, 


is determined using (equation 8) defined above. B;; values are determined experimentally from the particle size 
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analysis of the product of short-term grinding of a single starting mass of size j (BIJ method) described by [45], 


using (equation 12): 


1-P(0) 
| 
eee (12) 


B,. 
: 1— P,,,(0) 
i eee 
1-P.,.(f) 


P(t) is the weight fraction below the top size of size range i after grinding time t. 


3.2.3. Back calculation parameters 


The values of the optimized parameters of the C, S and B matrix are determined using the back-calculation 
procedure described in Fig. 3 below. Back-calculation, a technique which involves calculating the model 
parameters that best fit the experimental data, is also widely used and offers significant advantages over direct 
measurements [34],[47],[53]. Linear and non-linear optimization techniques are usually combined, or combined 


with other approximate solutions, to recalculate the model parameters [54],[55]. 


Figure 3 shows a diagram of the procedure for retro-calculating the parameters of the selection, grinding and 
classification functions. 


Numerical 
Data 


Matrix PBM SSR / \»|Back - calculated 
parameters equation re , / parameters 
Experimental rf 
Data 


Figure 3. Schematic diagram of the back — calculation procedure adapted from [24] 


Experimental 
Data 


The optimal values of the fitted parameters are estimated by minimizing the objective function (sum of squares of 


the residuals), defined by (equation 13) below: 


Min SSR = yy, 


i=l 


Ying) (13) 


xp, 1 = 


The optimization method to back-calculate these parameters of matrix S, B and C was carried out using the 
"fmincon" nonlinear optimizer in MatLab'™ 2019. "fmincon" uses the "point-interior" algorithm to solve the 


nonlinear optimization problem in question, and is suitable for general nonlinear optimization [47],[56]. 
3.2.4. Specific energy consumption modelling 


For our study and based on the work of [5],[21],[59], we hypothesize that the specific energy consumed E,,.. by the 
local fabricated hammer mill is modeled bt the Rittinger energy law, whose expression is described by the equation 


(14): 


1 1 
Ek 14 
C 7q> (14) 


spec 
2 1 
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with d, and d, are the feed and product size, respectively; and K is a Rittinger constant. In this equation (14), the 


product size d, is calculated by the equation (15): 


d, =) Xp; (15) 

i=l 
where pi is the mass fraction of product retained of sieve of size d,;, during the particle size analysis, and Xpi is the 
geometric mean diameter of particles of i” sieve, assumed to be described by the equation (16) below: 


xX .=,/D,[D 


pi pti p.i-l 


(16) 


where d,; and d, ;,; the diameter of sieve openings respectively of the i™ sieve and (i+1)* sieve. 


The value of the constant K is determined by adjusting the values of specific energy consumed obtained 


experimentally with those simulated using the relationship (14) of the rittinger’s energetic law. 
3.2.5. Determining optimum operating conditions 


In our study, we wish to have only particles whose size is greater than 0.25 mm; the standard stipulates that less than 
25% of the flour must pass through a sieve with a mesh size of 0.25 mm [57] so our study aims to minimize particles 


with sizes smaller than 0.25 mm when grinding cereal corn in a locally manufactured hammer mill. 


In the case of a multi-criteria problem, the challenge is to determine a solution that is an acceptable compromise 
between all the criteria considered. In the case of grinding, we have defined two optimization criteria: the specific 
energy consumed and the quality of the ground product. Cereal grinding using a hammer mill generates large 
quantities of fine particles (less than 0.25 mm in size), which are harmful to the environment and the surrounding 
population. The quantity of fine particles generated depends on several parameters, including rotation speed, feed 
rate, sieve size and the characteristics of the material to be ground. Our study will focus solely on the "sieve size" 
variable, as its impact on specific energy and flour quality (particle size distribution and geometric mean size) has 


been shown to be predominant [9],[58]. 
Then, the functions to be optimized are as follows: 


n-2 


min p, = P,1+P, =1->P, (17) 
i=] 
min E ee (18) 
spec r d, d, 


The particle size distribution of outlet p is described by (equation 4) previous with the fragmentation matrix X 


described in (equation 5) previous. 


with constraints: 
pe Viele syn 19) 
P=! (20) 


Px = Pr + ?P, < 0,2 (21) 
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be optimized in our problem are the speed of rotation and the size of the sieve inserted in the mill. We note [Xopt] 


the optimal values of our parameters. 
== 4, Results and Discussion 
4.1. Experimental values of matrix S and B 


The experimental values of the Si selection matrix corresponding to maize particle sizes for different classes are 


shown in fig. 4. 


This shows the influence of the rotation speed on the milling frequency of maize grains. The selection function S 
increases with the rotation speed of the rotor. We also notice that when the feed rate is growth, the values of Si 
increase too. This suggest that a higher fragmentation of the particles needs a higher rotation speed with an 
important feed rate at the entrance. The theory of the particles rupture shows that the bigger particles break more 


easily because they contain microcracking more important than the smaller ones [60]. 


0.0120 
0.0115 
—_ — —— 2 
o ee ee 
__ 0,0110 eo 
= 
‘” 
~— 00105 + 
wo 
u— 0,0100 4 pp 
5 r 
@ 
0.0085 + 
= -¥ 
rf * 
> o,0000 4 -— 
_~ —e— TO 2000 term 
< TO & 3000 term 
0.0085 4 v ——-¢—- PO&2000 tim 
—i— PO 43000 trimn 
0.0089 +, -_ er 
4 6 8 10 12 14 16 


Size particles (mm) 
Figure 4. Influence of the feed rate and the rotation speed on the experimental values of S; 
The values of matrix S parameter obtained using experimental data is presented in table 4 below: 


Table 4. Parameters values of function S 


Opening hopper | Parameter a he 
Value 3.16 2.375 
PO 
Sdt error 0.068 0.0529 
Value 3.241 2.304 
TO 
Sdt error 0.056 0.0540 
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Experimental values for the grinding matrix are calculated using the BIT method (equation 12), which is used for a 


very short grinding time based on data obtained in experiment A. The results obtained are presented in Tables 5 and 


6 for the partially open and fully open hopper experiments. 


Table 5. Experimental values of Bj; and b,j at PO hooper 


V = 2000 tr/mn V = 3000 tr/mn 
Class of 
particles I IT Tl I IT Til 
aa Bi bit Bi2 bi2 Bi3 bi3 Bi bi1 Biz bi2 Bi3 bi3 
1 1 0 1 0 
2 1 0.216 1 0 1 0.159 1 0) 


3 | 0.784/ 0.239} 1 |0.357} 1 O |0.841/0.217) 1 {0.247}; 1 0 


4 | 0.545 | 0.082} 0.643}0.156} 1 | 0.451] 0.624/ 0.107] 0.753/0.158} 1 | 0.366 


5 | 0.463) 0.211) 0.487} 0.18 | 0.549} 0.178} 0.517] 0.144] 0.595 | 0.21 | 0.634 | 0.196 


6 | 0.252} 0.058 | 0.307 | 0.061 | 0.371 | 0.098 | 0.373 | 0.116 | 0.385 | 0.193 | 0.438] 0.151 


7 | 0.194) 0.112} 0.246} 0.154) 0.273 | 0.127} 0.257] 0.114] 0.192} 0.09 | 0.287 | 0.142 


8 | 0.082 | 0.056} 0.092 | 0.058 | 0.146 | 0.091 | 0.143} 0.046} 0.102 | 0.018 | 0.145} 0.05 


9 | 0.026} 0.017} 0.034} 0.018 | 0.055 | 0.034 | 0.097 | 0.082] 0.084 | 0.048 | 0.095 | 0.05 


10 | 0.009 | 0.009 | 0.016 | 0.016 | 0.021 | 0.021 | 0.015 | 0.015 | 0.036 | 0.036] 0.045 | 0.045 


The grinding matrix represents the cumulative size distribution of the child particle fractions resulting from the 


primary grinding of the mother particle. 


Table 6. Experimental values of b;; and b;; at TO hooper 


V = 2000 tr/mn V = 3000 tr/mn 
Classes of 
particles I II III I II Il 
sizes 
Bi bit Biz biz Bi3 bi3 Bi bit Biz biz Bi3 bi3 
1 | 1.000] 0.000 1.000 | 0.000 
2 | 1.000} 0.205} 1.000} 0.000 1.000} 0.109} 1.000} 0.000 
i= 
3 | 0.795] 0.233} 1.000} 0.160} 1.000} 0.000} 0.891] 0.207] 1.000] 0.247] 1.000] 0.000 
4 | 0.562] 0.278 | 0.840} 0.190} 1.000} 0.170} 0.684] 0.397] 0.753] 0.244] 1.000] 0.284 
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5 | 0.284} 0.111] 0.650 | 0.274 | 0.830) 0.176 | 0.287) 0.114} 0.510} 0.202) 0.716) 0.179 


6 | 0.173} 0.088 | 0.376 | 0.225 | 0.654 | 0.256 | 0.173 | 0.088 | 0.308 | 0.139 | 0.538 | 0.297 


7 | 0.085 | 0.032 | 0.151 | 0.055 | 0.398 | 0.227 | 0.086 | 0.032 | 0.169 | 0.086] 0.241] 0.135 


8 | 0.053] 0.022 | 0.096 | 0.044 | 0.171 | 0.097 | 0.054 | 0.009 | 0.083 | 0.040 | 0.106] 0.005 


9 | 0.031} 0.016 | 0.052 | 0.031 | 0.074 | 0.041 | 0.045 | 0.029 | 0.057 | 0.031 | 0.101 | 0.017 


10 | 0.015] 0.015] 0.021] 0.021) 0.033] 0.033 | 0.016] 0.016] 0.026) 0.026) 0.084) 0.084 


The values of B;; as a function of particle size obtained from BII calculations for the parent fraction > 10 mm in size 


are shown in Fig. 5 above. 
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Figure 5. Experimental values of B;; for maize size fraction > 10 mm 


The following observations can be made from Fig. 5 above, when compared with the results obtained from other 


maize grain fractions: 


e The feed rate and rotation speed slightly improve the size distribution of the ground product; finer particle sizes 
are observed when the hammer mill is sufficiently fed and rotates faster, showing that these two conditions favor a 
greater quantity of grain being subjected to fractionation; and the rotation speed being high, the fractionation energy 


is considerable, resulting in smaller particle sizes. 


e The results obtained show that changes in mill operating conditions do not lead to significant changes in the 
values of the parameters @ and y. This observation is in agreement with previous studies [61],[62] which also 


indicate that the parameters b;; are independent of mill operating conditions. 
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These plots show that the grinding function does not depend on the initial particle size, so the grinding distribution 


function can be considered normalizable. In this way, the values of B can be fitted to the empirical relationship, and 


the parameters of the grinding function are determined and presented in Table 7 below: 


Table 7. Parameters values of function B 


Opening 
hopper Parameters ct) B y 
Value 0.328 | 3.102 | 0.809 
PO 
Sdt error 0.004 | 0.0141 | 0.008 
Value 0.331 | 3.275 | 0.812 
TO 
Sdt error 0.004 | 0.0129 | 0.007 


4.2. Result of experiment B and estimated parameter values of the classification function 


The same experiments were then carried out, but the sieve was left open. The values obtained after particle size 
analysis of the flours obtained are given in the table below. The cumulative mass fractions obtained from 


experiment B are presented in Table 8 below for grain sizes ranging from 6.3 to 10 mm. 


Table 8. Cumulated fractions obtained during experience B for the partially opened hopper 


V= 2000 tr/mn V= 2000 tr/mn V= 3000 tr/mn V= 3000 tr/mn 


ds = 2 mm ds = 8mm ds = 2mm ds = 8mm 


99.52] 98.92 ie 95.57] 96.37 | 95.75| 99.31 97.89 


a a a a a 


This table 8 shows the influence of rotation speed and initial particle size on the particle size distribution at the 
outlet: for a rotation speed of 2000 tr/mn, the dso (the mean or average particle size of particles) of the products is 
less than 0.63 mm when using a 2 mm sieve it also mean that 50% of the particles are smaller than 0.63 mm, but is 
greater when using an 8 mm sieve. At speeds of 3000 tr/mn, the dso is less than 0.63 mm. This shows that as hammer 


speeds increase, the average geometric size of the final products decreases. 
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From the results obtained below and using the formulas defined in (equation 9) and (equation 10), we obtain the 


parameter values shown in Table 9. 


Table 9. Parameters values of the classification matrix 


Parameters Oy Oy 03 4 


Value 1.097 | 0.844 2.133 1.243 


Hopper 
PO Stddev. | 0.088 | 0.0704 | 0.0618 | 0.0413 
: Value 1.123 | 0.873 | 2.152 | 1.370 
opper 
TO 


Std dev. 0.092 | 0.0711 | 0.0631 | 0.0441 


4.3. Matrix parameters values after back - calculation 


The parameters of the model (A, 42, 9, y, B, @1, @2, @3 and a4) vary as a function of the process variables. Therefore, 
the idea here, is to find a single set of parameters describing all the experimental data consisted in performing 
different preliminary estimations and to correlate their parameter values obtained to the process conditions. To 
increase accuracy, the initial parameters of the selection, breakage distributions and classification functions were 


initially fixed turn-by-turn to reduce the time of the back-calculation. 


Table 10. Optimized function parameter values 


Opening 
Parameter my ho 0) B Y Oly On 03 Oly 
hopper 
Value 3.05 | 2.30 | 0.330) 3.28 | 0.811) 1.095 | 0.866 | 2.141 | 1.31 
PO 
Std dev 0.061 | 0.055 | 0.004 | 0.012] 0.008] 0.089 | 0.065 | 0.062 | 0.047 
Value 3.12.) 2.28 |0:330)| 3.35 0811) 1-105) 0.859 2.146 | 142 
TO 
Std dev 0.064 | 0.057 | 0.004 | 0.014} 0.008] 0.083 | 0.070 | 0.062 | 0.046 


Table 10 summarizes the best combination of back-calculated parameters obtained after simulation with data from 


Experiment C. 


Fig. 6 shows the model's predictive behavior for corn grinding for trials 1, 2, 3 and 5. It compares the experimental 
size distributions obtained in Experiment C with the simulated distributions for the identified operating cases. In 
this figure, we observe very good agreement between experimental and predicted values, within the limits of 


experimental error, as confirmed by the value of the coefficient of determination R’. 
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Figure 6. Comparison between experimental and simulated size distribution for partially opened feed hopper (PO) 


Overall, Fig. 6 shows a good correlation between the model and experimental data. However, there is a slight lack 
of precision in the case of higher rotational speed and larger sieve aperture size. In the most favorable cases, the 
quality of the prediction is visually excellent, to the extent that the prediction error is of the same order of magnitude 


as the experimental measurement errors. 
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Figure 7. Comparison of experimental and simulated size distribution for a partially open feed hopper (PO) 
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Experimental and simulated particle size distributions for trials 7, 8, 9 and 10, corresponding to higher mill feed 


rates (hopper is fully open), are shown in Fig. 7, together with the coefficient of determination (R’). 


These graphs show strong agreement with the experimental data, particularly for test 9 performed at a rotational 
speed of 2000 rpm with a sieve opening size of 8 mm; and for test 8 performed at a rotational speed of 2500 rpm 


with a sieve opening size of 5 mm. 


For trial 10 (3000 rpm with 8 mm sieve opening size), the correspondence is weaker, although these are the variable 
parameters for which less experimental data were available; this observation was also made when the mill feed rate 
was lower, although these were the variable parameters for which less experimental data were available ; this 
observation was also made when the mill feed rate was lower, suggesting that the model would not fit when for a 
very high machine feed rate a large screen with a higher rotor speed is used. In general, however, as shown in Fig. 7, 


the model is able to correctly predict changes in sieve opening size, as well as changes in rotational speed. 


From this study, it can be noted that the grinding matrix was independent of rotor speed, mill feed rate and sieve 
size, but was dependent on grain characteristics. The selection and classification of the material inside the mill 
depended on the physical configuration and operating characteristics of the mill, as well as on the size of the input 


particles. 
To improve the predictive quality of our model, the following observations are made: 


- Itseems very important to feed the model with an appropriate set of data for its calibration; this can be achieved 


by increasing the number of experimental points. 


- Itis also likely that the nature of the material affects the quality of the prediction. The model specifically models 
the grinding of maize kernels of a particular variety with a fixed moisture content; characterization errors may 
therefore be found to be caused by the relative humidity of the working environment, which creates or simulates 


agglomerations of fine particles. 


Despite the good correlation between experimental and simulated observations or data, it is possible that a more 
mechanistic model could improve the applicability of the study [9],[63]. However, it should be noted that particle 
crushing or fragmentation is the result of a number of microphenomena that are not very well understood. Complex 
interactions between particles, particle stress conditions and the environment will determine the outcome of these 


industrial processes. 
4.4. Optimum operating conditions for a locally manufactured hammer mill 


For a better presentation of the optimization problem, we assume that the particle size distribution of the output 


product as a function of the input product is defined by equation (22). Equation (22) is deducted from (4). 
p=[G]m (22) 


where G is a lower triangular matrix, whose elements is defined in equation (23): 


[G, Joc =f (N.ds) (23) 


izi 
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n-l n 
max [2 ->)°G,_,,(N,ds\im, - 1G, (N, asym 
i=l 


i=l 


(24) 
1 1 


Spit: G,lim, “ 
j=l i=l 


min kK 


Table 11 below shows the results obtained on Matlab using the "gamultiobj" multi-objective optimization 


algorithm. 


Table 11. Optimum solutions for [Xopt] values in the TO hopper 


Parameters N ds 
Values 2750 tr/min 4mm 
Std dev 25 0.97 


From the values obtained in Table 11, we simulated the particle size distribution of the product and calculated the 


specific energy consumed. The results are shown in Table 12 below: 


Table 12. Simulated flour mass fractions and specific energy consumption for optimum mill operating conditions 


Espec Size of Sieve (mm) 
dso (mm) 
(KWh/t) 1.6 | 1.00} 0.8 | 0.63 | 0.425| 0.25) pan 
[Xopt] 4.48 0.12] 0.13] 0.15) 0.19} 0.23 | 0.13] 0.04 763.06 
1,0 
a. 
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Figure 8. Simulated particle size distributions for [Xopt] compared with tests 5 and 10 
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Calculating the mass fraction of the last two analysis sieves and comparing with those obtained from experiment C, 


we note a 20% reduction in fine particles, and 5% lower energy consumption. 


Calculating the mass fraction of the last two analysis sieves and comparing with those obtained from experiment C, 


we note a 20% reduction in fine particles, and 5% lower energy consumption. 


Fig.8 above shows the simulated particle size distribution for the optimum operating conditions, which significantly 
reduce fine particle production with low energy consumption. It can be seen that, in the case of our mill, it would be 
preferable to operate the mill at a higher feed rate, as this would help to reduce the production of fine particles, but 


with a sieve with aperture sizes close to 4 mm. 


Although there is a difference between partially open and fully open hopper feeding, this difference is not 


significant in terms of values. 
= 5. Conclusion 


Typically, approaches to milling optimization have focused more on improving milling efficiency in terms of 
quantity produced and energy consumed, thus aiming for higher quantity produced and lower energy consumption 
during the process. This left very little room for the use of modelling to reduce the quantity of fine particles which 
are generally lost during transport of the flour produced and also constitute a health risk for the operators of this 
equipment. The aim of this work is therefore to use milling modeling to optimize flour size distribution, so that the 
quantity of fine particles produced during milling is reduced or minimized, without increasing the system's energy 
consumption. To this end, the population balance approach was used to model the particle size distribution of the 
flour produced; the model obtained proved satisfactory for predicting the particle size distribution of the products 


leaving our locally-produced hammer mill. 


In this study, we have shown that the selection function for particle fracture and the classification function are 
closely related to rotation speed and sieve size. Both depend on the speed at which particles are available for impact 
fracture or to pass through the sieve aperture size. As rotational speed increases, more particles can be selected for 
fracture and pass through the sieve, reducing mass retention. Overall, a model representing the grinding operation in 
hammer mills has been developed and validated on the basis of steady-state experimental data from mills of 
different scales. These results are of great value, as they indicate that it is possible to use the fitted model to 


confidently predict new operating points. 


Next, we studied the operation of our mill by defining the optimum operating points of our mill; it emerged that the 
optimum rotational speed Nop = 2750 rpm combined with a sieve size don = 4 mm, contribute to producing fewer 


fine particles with low energy consumption. 


Looking ahead, we can carry out further tests to validate our model on other types of hammer mills. In addition, 
other more mechanistic functions can also be exploited to broaden the applicability of the proposed model if they 
can take into account, in addition to material properties, particle size and process energy consumption. This will 


also enable the model to be more accurate beyond the limits of experimental conditions. 
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